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ABSTRACT 

The first challenge in the formation of both terrestrial planets and the cores 
of gas giants is the retention of grains in protoplanetary disks. In most regions 
of these disks, gas attains sub-Keplerian speeds as a consequence of a negative 
pressure gradient. Hydrodynamic drag leads to orbital decay and depletion of 
the solid material in the disk, with characteristic timescales as short as only a 
few hundred years for meter-sized objects at 1 AU. In this paper, we suggest 
a particle retention mechanism which promotes the accumulation of grains and 
the formation of planetesimals near the water sublimation front or "snow line." 
This model is based on the assumption that, in the regions most interesting for 
planet formation, the viscous evolution of the disk is due to turbulence driven 
by the magneto-rotational instability (MRI) in the surface layers of the disk. 
The depth to which MRI effectively generates turbulence is a strong function of 
the grain size and abundance. A sharp increase in the grain-to-gas density ratio 
across the snow line reduces the column depth of the active layer. As the disk 
evolves towards a quasi-steady-state, this change in the active layer creates a local 
maximum in radial distribution of the gas surface density and pressure, causing 
the gas to rotate at super-Keplerian speed and halting the inward migration of 
grains. This senario presents a robust process for grain retention which may aid 
in the formation of proto-gas-giant cores preferentially near the snow line. 

Subject headings: planetary systems: protoplanetary disks - solar system: for- 
mation - turbulence 



1. Introduction 

More than 200 extrasolar planets have been discovered with radial velocity surveys, 
corresponding to 10 to 15% of the sample of mostly nearby solar-type stars(Marcy et al. 
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2005). In the current paradigm, formation of both terrestrial and gas giant planets begin 
with the collisional growth of planetesimals, the latter only accreting gas once the solid 
cores have reached several earth masses (Pollack et al. 1996). The growth rate of these 
protoplanetary embryos is determined by the surface density distribution of planetesimals, 
S p , in their nascent disks. In the typical construction of planet-formation sequence (eg. 
Kokubo & Ida 2002, Ida & Lin 2004, hereafter IL04), authors have often utilized an empirical 
minimum mass nebula model" (MMNM) in which E p ~ 10r A u /2 g cm , where r^u is the 
radial distance from the central star. This ansatz is based on the assumption that the planets 
in the solar system locally retained all the heavy elements and a fraction of the gas in the 
primordial solar nebula (Hayashi et al. 1985). Within the observationally inferred range of 
protostellar disk masses (Beckwith & Sargent 1996), gas giants can readily form near the 
snow line on disks with 2 to 10 times enhancement over MMNM(IL04). 

Despite the successful reproduction of the planets' observed mass-period distribution, 
the global scaling of E p with that of the MMNM is questionable. The main uncertainties 
are: (1) whether heavy elements can be efficiently incorporated into planetesimals which 
do not migrate significantly and (2) whether there is a sufficient E p enhancement near the 
snow-line to favor the emergence of gas giant planets. The observed surface density of the 
dust, Ed, is in the form of grains ranging from a few /iin to mm. These entities must grow 
into embryos larger than 10 3 km in gaseous disks. In most regions of the disk, the pressure, 
P, decreases with r so that the azimuthal velocity of the gas, V^, is sub-Keplerian. Grains 
larger than a few cm are decoupled from the gas and move at Keplarian speeds so they 
experience head winds and undergo orbital decay on time scales smaller than the growth 
time scale by a factor Z = E^/E 9 (Garaud et al. 2004). In principle, grains at a few AU can 
be retained and then converted into planetesimals either through sedimentation followed by 
gravitational instability (Goldreich & Ward 1973; Youdin & Shu 2002; Garaud & Lin 2004) 
or through cohesive collisions (Supulver & Lin 2000) after the photoevaporation of disk gas 
increases Z above unity. However, the severe depletion of the disk gas will also limit the 
supply for the gas giants to accrete their massive envelopes. 

In this letter, we propose a grain retention mechanism which can promote rapid planet 
formation around sublimation fronts in protoplanetary disks. We assume the disk evolution 
during the classical T Tauri phase to be regulated by a magneto-rotational instability (MRI; 
Balbus & Hawley 1991). However, MRI only operates in regions with sufficient ionization 
fraction to couple the gas to the magnetic field. The inner regions of typical protostellar disk 
are thermally ionized (Pneuman and Mitchell 1965; Umebayashi 1983). At greater distances 
from the disk center (r > 0.1 AU), stellar x-rays and diffuse cosmic rays ionize the surface 
layers of the disk down to a certain column density, E^ (Glassgold et al. 1997). Thus, the 
disk may have a layered structure in which the MRI-induced accretion flow occurs on its 
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active surface which surrounds a highly neutral and inactive "dead zone" (Gammie 1996). 
The magnitude of depends strongly on the relative abundance and size distribution of 
grains (Wardle & Ng 1999; Sano et al. 2000). Across the snow line or the dust destruction 
front, there are significant changes in the amount of solid material in the form of grains. In 
a disk with a solar composition, the total mass fraction in grains decreases by at least 50% 
due to the sublimation of water ice interior to the snow line (Lodders 2003). Efficient inward 
transport of solid materials combined with the slow diffusion of the heavy elemental vapor 
can further enhance the contrast in grains' surface density (Cuzzi & Zahnle 2004). This rapid 
radial change in dust properties and the corresponding rapid change in the column density 
in the active region leads to local maxima in S 9 and P in a steady-state disk. As noted 
by other authors (Bryden et al. 1999; Haghighipour & Boss 2003), this pressure-gradient 
reversal induces a super-Keplerian V$ and causes solid material to build up near the local 
pressure maxima. 

In order to highlight the dominant process we adopt a simple, 1-d, steady-state model 
for the gas distribution, in which the gas mass flux in the radial direction, M g = 2^H g V r r 
(where V r is the radial velocity of the gas), is independent of r. Interior to 5 AU in a typical 
protoplanetary disk the disk should rapidly relax to this quasi-steady-state (Beckwith & 
Sargent 1996). In §2, we present a simple analytical model for the structure of a disk 
with variable viscosity and describe the conditions under which these changes may induce 
a pressure-gradient reversal. In §3, we present numerical results to show how these changes 
manifest themselves around the dust destruction front and the snow line, at different stages 
in the disk evolution. Finally, we summarize our results and discuss their implications in §4. 



2. Model Description 

We will parametrize our model in the framework of the familiar ad hoc a-prescription in 
which the efficiency of angular momentum transport is approximated by an effective viscosity 
v = a(r)c s h, where c s and h = c s /Qk refer to the sound speed at the mid plane and the 
isothermal density scale height, respectively (Shakura & Sunyaev 1973). However, we deviate 
from the standard model in that we allow a to vary with r. The Keplerian frequency around 
a star with mass M* is VLk = (GM^/r 3 ) 1 ^ 2 . The surface density at each radius in a disk with 
mass accretion rate M g is £ 9 = M g /(37ru(r)) (Pringle 1981). 

The interaction of the gas with the sedimented grains (super-cm size) leads to their 
orbital evolution. The speed of their migration depends on their sizes, and the radial velocity 
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of the maximally affected particles is 



2pdr K ' 



(see Weidenschilling (1977) for derivation). In most regions of the disk P monotonically 
decreases with r so that < VLkt and gas drag causes the grains to undergo orbital 
decay. However, if dP/dr > then v r > 0, and particles will drift outwards rather than 
inwards. If we assume a vertically isothermal disk, the pressure at the midplane is given by 
P = Y, g cj./(27rh). If we assume a power-law temperature distribution in the radial direction, 
T = T r^y, and use the steady-state S 9 , then v r from eq(0Q) will be positive if 
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We can approximate a rapid radial variation in the viscosity at radius tq as a smoothed step 
function in a, 
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where Ar is the characteristic width of the transition zone. For this simple a(r) prescription, 
the particle outflow condition is 



1 



Ar / 1 Act 
V < \ln a (3-q/2)' 



(4) 



where Aa = a.\ — a>2 and a = (cti + a2)/2. 



In a layer-accreting disk there is physical motivation to describe the variations in a with 
respect to the depth of the region unstable to MRI turbulence. In this case a is the vertically 
integrated value including both the efficient angular momentum transport from the active 
layers in the disk and a small residual contribution ao from the MRI stable region such that 



a 



-ttMRI + a , 



(5) 



where cimri = 1-8 x 10~ 2 (/9/1000) _1 is the appropriate scaling for MRI turbulence (Sano et 
al. 1998). For the plasma parameter (/5 = P/P mag ) we use (3 = 1000. For disks with weak 
magnetic fiends, reasonable values for (3 are on the order of 100 to 1000 (Sano et al. 2000). 

Although an actual detailed description of is not important for the effect we consider 
here, it represents some mechanism of angular momentum transport in the "quiescent" part 
of the disk beyond the dust-destruction front (e.g. convective instability in the vertical 
direction (Lin & Papaloizaou 1980), gravitational torques (Laughlin & Rozycza 1996), linear 
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Rossby wave instability (Li et al. 2000), baroclinic instability (Klahr & Bodenheimer 2003; 
Klahr 2004), linear stratorotational instability (Dubrulle et al. 2005; Shalybkov & Riidiger 
2005), etc.) so that it is possible for the disk gas to evolve into a quasi-steady-state, in which 
M g becomes an essentially independent function of r. In the absence of any other angular 
momentum transport mechanism (i.e. «o = 0), a steady state is unattainable and gas will 
continue to accumulate until another instability is generated. A pressure-gradient inversion 
can still develop, albeit it will be non-trivial to retain the accumulated grains if the disk 
becomes gravitationally unstable. 



The steady-state surface density profile is: 

M 



«MRI 



(6) 

The column depth of the region unstable to MRI turbulence can locally be described as a 
powerlaw XU oc r p when the total column depth S > E^. At a sublimation front in the 
layer-accreting region of the disk, the active region may be described locally as 



M,2; 
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and XU,2 are the column depth at that radius assuming the the two different particle 
populations on either side of the transition. In this case the particle outflow condition is 
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where Emri = M g / (37raMRic s /i) is the steady-state surface density if the entire disk were 
MRI active. In the case when Smri 3> this solution reduces to 



Ar 
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(9) 



r V 7T E MRI (3 - q/2) ' 

Note that this is independent of ao, the magnitude of the residual viscosity. However, it is 
only meaningful in situations where a is large enough that the disk can reach a steady-state 
solution without becoming unstable. 



3. Results 

3.1. Dust Destruction Front 

The modified step function in eq© can be most easily applied in the inner region of the 
disk near the dust destruction front. The gas is thermally ionized until T ~ 1000K, so the 
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dust-destruction front is roughly coincident with the outer boundary of the thermally ionized 
region. At this location, the disk makes a transition from being completely MRI-active to 
having a very large gas column density which is stable to MRI. As the active region switches 
from covering the entire thickness of the disk to being virtually non-existent in the presence 
of a significant population of small particles (Turner et al. 2007), we can utilize our simple 
a description from eq(j3j) where «i = «mri and «2 = «o- 

Due to the expected large decrease in viscosity, the condition stated in eqfll]) is easily 
satisfied at the dust-destruction front for even rather broad transition regions. Therefore, 
we use this clean case to demonstrate the manner in which the pressure gradient is affected 
by the viscosity in a steady-state disk. Fig. [1] shows the effect for a particular a profile with 
cto = 10~ 3 for different power-law temperature profiles. While the surface density changes 
according to the temperature profile, the pressure gradient is very weakly dependent on 
temperature. We also show the evolution of an initially well-mixed population of maximally 
migrating particles. While a realistic population of particles will not accumulate as quickly as 
this maximally drifting population, this is illustrative of the efficiency of particle retention at 
the pressure maxima. It is also important to note that the location of the pressure maximum 
is exterior to the sublimation front, meaning that the material trapped in this region will 
not be lost due to sublimation. 



3.2. Snow Line 

A reversal of the pressure gradient around the snow line will only occur at certain stages 
in the disk evolution. In an MRI-driven disk the structure of the disk at the snow line can 
be described in three stages. Initially, the column density at the snow line is very large, 
corresponding to a high steady-state mass accretion rate. The MRI-active layer is negligible 
due to the presence of numerous small grains. If the mass accretion rate is high (as in disks 
around FU Ori systems), £ 9 will increase to such large values that the disk will be unstable 
to gravitational instabilities (Laughlin & Bodenheimer 1994). The small contribution of 
the change in the active layer at the snow line will not be important in these early disks. 
Later, the mass accretion-rate decreases and/or the grains grow so that the structure of 
the MRI-active layer strongly affects the total disk column depth. In this stage there is 
a dip in the steady-state surface density profile corresponding to the increase in viscosity 
interior to the snow line. In these disks the pressure gradient is inverted, and solid material 
accumulates. Finally, the mass accretion rate, and therefore surface density, will decline to 
the point that the entire disk is MRI-active, and this trapping mechanism will no longer 
be effective. Regardless of the grain properties this stage will have occured by the time the 
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mass accretion rate drops to the level observed among the weak-line T Tauri stars. 

To demonstrate these stages Fig. [2] shows the resulting steady-state surface density pro- 
files for a series of accretion rates and illustrative examples of a® (not intended to constrain 
the possible range of «o). We approximate the size of the active layer from chemical equi- 
librium calculations to determine the linear stability to MRI turbulence (Ilgner & Nelson 
2006 (model 4); Neal Turner private communication). In these calculations we assume a 
grain size of 1/im. If the grain size is substantially smaller than this the column depth of the 
active region at the snow line may be too small to generate turbulence. The most obvious 
feature in the surface density profiles is the sharp decrease in the surface density at the 
edge of the dead- zone (beyond which S 9 = £4). The effect that we are interested in is the 
small change in surface density at the snow line, around 3 AU. The curves span the range of 
phases described abover, so only the middle two have pressure inversions at the snow line. 
We note that the assumed temperature profile with q = 0.5 is not self-consistent with the 
disk structure assuming a passively irradiated disk, but as the pressure gradient is only very 
weakly dependent on the temperature we neglect the self-consistent calculation at this time. 



4. Discussion and Conclusions 

We have demonstrated that, in a disk whose evolution is dominated by the MRI-driven 
turbulence, the sensitivity of the instability to the grain properties creates regions with radial- 
pressure-gradient inversion and results in the accumulation of migrating particles. This 
mechanism provides a potential resolution to the problem of grain retention in protoplanetary 
disks. These radial-pressure-gradient inversions occur at the sublimation fronts, such as the 
snow line, with a modest grain depletion as long as the transition occurs over a short radial 
extent. 

In the presence of another viscous mechanism in the MRI-stable zones the disk may be 
able to reach a quasi-steady-state. Under these conditions, the solid material will continue 
to flow into the pressure maximum until the amount of solid material at this radius becomes 
comparable to the amount of gas in the disk. Thereafter, the dust transfers angular mo- 
mentum to the gas in the inactive shielded layers, forcing it to rotate at the local Keplerian 
speed and to diffuse to larger radii. The region of dust accumulation then expands slowly 
outward. In the absence of another source of viscosity, a steady-state may not be attainable. 
In this case, gas will continue to accumulate, maintaining a pressure maximum which will 
trap solid material until other dust or gas instabilities start to dominate the disk evolution. 

The particle collected at the snow line should rapidly grow through cohesive collisions. 
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At the snow line, the grains of all sizes are likely to be covered with marginally molten 
frost which can promote the cohesive nature of the collisions (Bridges et al. 1996). The 
particles may also be able to grow past the "meter- sized barrier" (Weidenschilling 1984; 
Cuzzi & Weidenschilling 2006). First, as at the pressure maximum all material is moving 
at Keplarian velocity, large objects will no be "sand-blasted" by smaller grains. If particles 
grow past the meter-sizescale to the regime where they are not well coupled to the turbulent 
flow there is no further barrier to their growth as there would be at other poisitions in the 
disk. Additionally, as the particles are collected in regions where the midplane is stable to 
MRI-turbulence, the question of whether meter-sized objects (coupled to the largest eddies) 
will undergo destructive collisions depends on the properties of whatever mechanism for 
angular momentum transport may occur in the "dead zone." If the disk is weakly turbulent 
(cto < 10~ 4 ) or if the trasport of angular momentum is due to a mechansism not involving 
turbulence then even this aspect of the "barrier" can be overcome. 

We thank Neal Turner for many enlightening discussions and for providing us with the 
calculations of the depth of the MRI-active layer. We also thank our referee, John Cham- 
bers, for his helpful comments. This work has been supported by NASA (NNG05G-142G, 
NNG04G-191G, NNG06-GH45G), JPL (1270927), NSF(AST-0507424) and IGPP1326. 
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Fig. 1. — The effect of an example variable a(r) shown in panel (a) on the gas properties. 
Panel (b) shows the corresponding steady-state gas distribution with temperature profile T = 
T r^fj. The curves represent g=0.25, 0.5, and 0.75 (dotted, solid and dashed, respectively). 
Panel (c) shows the corresponding pressure gradient. Panel (d) shows the evolution of an 
initially well-mixed solid population after 500 orbits at 0.1 AU (grey lines show every 100 
orbits). 
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Fig. 2. — The steady-state profiles for two different alphas and different mass accretion rates 
M s =10~ 9 , 3 x 10" 9 , 8 x 1(T 9 , 3 x 1(T 8 A^ear" 1 (black dot-dashed, dashed, solid and dotted, 
respectively). The grey dashed line shows the Toomre stability criteria. The grey solid line 
shows the with \[im grains with a dust-to-gas ratio changing from 0.005 to 0.011 across 
the snow line (Lodders 2003). Panels (i)-(iv) show the corresponding pressure gradients for 
each M g , which are independent of a . 



